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A differential Lyapunov framework 
for contraction analysis 

F. Forni, R. Sepulchre 
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^jq[ Abstract 

Lyapunov's second theorem is an essential tool for stability analysis of differential equations. The 
paper provides an analog theorem for incremental stability analysis by lifting the Lyapunov function to 
the tangent bundle. The Lyapunov function endows the state-space with a Finsler structure. Incremental 
1 stability is inferred from infinitesimal contraction of the Finsler metrics through integration along 

solutions curves. 

1/3 '■ 

l ; I. Introduction 

>: 

At the core of Lyapunov stability theory is the realization that a pointwise geometric condition 

"\t" ■ 

Q\ 1 is sufficient to quantify how solutions of a differential equation approach a specific solution. The 

<n: 

geometric condition checks that the Lyapunov function, a certain distance from a given point 
to the target solution, is doomed to decay along the solution stemming from that point. By 
integration, the pointwise decay of the Lyapunov function forces the asymptotic convergence 

>■ 

to the target solution. The basic theorem of Lyapunov has led to many developments over 
the last century, that eventually make the body of textbooks on nonlinear systems theory and 
nonlinear control Ifl31 , [1311 , IfTOl , [fTTTl . Yet many questions of nonlinear systems theory call for 
an incremental version of the Lyapunov stability concept, in which the convergence to a specific 
target solution is replaced by the convergence or contraction between any pairs of solutions (3), 
lfT8l . Essentially, this stronger property means that solutions forget about their initial condition. 
Popular control applications include tracking and regulation ll24l . Il2ll . observer design , 
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[|28"1 , coordination, and synchronization [38], to cite a few. Those incremental stability questions 
are often reformulated as conventional stability questions for a suitable error system, the zero 
solution of the error system translating the convergence of two solutions to each other. This 
ad-hoc remedy may be successful in specific situations but it faces unpleasant obstacles that 
include both methodological issues - such as the issue of transforming a time-invariant problem 
into a time-variant one - and fundamental issues - such as the issue of defining a suitable 
error between trajectories -. Those obstacles also apply to the Lyapunov characterizations of 
incremental stability that have appeared in the recent years, primarily in the important work of 
Angeli 0- 

In a seminal paper ifTBl . Lohmiller and Slotine advocate a different angle of attack for nonlinear 
stability analysis. Their paper brings the attention of the control community to the basic fact 
that the distance measuring the convergence of two trajectories to each other needs not be 
constructed explicitly. Instead, it can be the integral of an infinitesimal measure of contraction. 
In other words, the often intractable construction of a distance needed for a global analysis can 
be substituted by a local construction. At a fundamental level, this approach brings differential 
geometry to the rescue of Lyapunov theory. The contraction concept of Lohmiller and Slotine 
- sometimes called "convergence" in reference to an earlier concept of Demidovich ][22| - has 
been successfully used in a number of applications in the recent years 0, ][24|. ETll . []38l . Yet, 
its connections to Lyapunov theory have been scarse, preventing a vast body of system theoretic 
tools to be exploited in the framework of contraction theory. 

The present paper aims at bridging Lyapunov theory and contraction theory by formulating 
a differential version of the fundamental second's Lyapunov theorem. Assuming that the state- 
space is a differentiable manifold, the classical concept of Lyapunov function in the (manifold) 
state-space is lifted to the tangent bundle. We call this lifted Lyapunov function a Finsler- 
Lyapunov function because it endows the differentiable manifold with a Finsler structure, which 
is precisely what is needed to associate by integration a global distance (or Lyapunov function) 
to the local construction. We formulate a Lyapunov theorem that provides a sufficient pointwise 
geometric condition to quantify incremental stability, that is, how solutions of differential equa- 
tions approach each other. The pointwise decay of the Finsler-Lyapunov function in the tangent 
space guarantees that a suitable (integrated) distance function decays along solutions, proving 
incremental stability. 
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There are a number of reasons that motivate the Finsler structure as the appropriate differential 
structure to study incremental stability. Primarily, it unifies the approach advocated by Slotine 
- which equips the state-space with a Riemannian structure - and alternative approaches to 
contraction, such as the recent approach by Sontag ll33l based on a matrix measure for the 
local measure of contraction. Examples in the paper further suggest that the Finsler framework 
will allow to unify the application of contraction to physical systems - typically akin to the 
Riemannian framework of classical mechanics - and to conic applications - typically akin to 
(non-Riemannian) Finsler metrics - such as consensus problems or monotone systems encoun- 
tered in biology. 

A primary motivation to study contraction in a (differential) Lyapunov framework is to make 
the whole body of Lyapunov theory available to contraction analysis. This is a vast program, 
only illustrated in the present paper by the very first extension of Lyapunov theorem based on 
LaSalle's invariance principle. Although we are not aware of a published invariance principle for 
contraction analysis, its formulation in the proposed differential framework is a straightforward 
extension of its classical formulation and we anticipate this mere extension to be as useful for 
incremental stability analysis as it is for classical Lyapunov stability analysis. 

We also include in this paper an extension of the basic theorem to the weaker notion of 
horizontal contraction. Horizontal contraction is weaker than contraction in that the pointwise 
decay of the Finsler-Lyapunov function is verified only in a subspace - called the horizontal 
subspace - of the tangent space. Disregarding contraction in specific directions is a convenient 
way to take into account symmetry directions along which no contraction is expected. This 
weaker notion of contraction is adapted to many physical systems and to many applications where 
contraction theory has proven useful, such as tracking, observer design, or synchronization. Those 
applications involve one or several copies of a given system and only the contraction between 
the copies and the system trajectories is of interest. 

The rest of the paper is organized as follows. The notation is summarized in Section lU 
Sections [Till 135] |V] contain the core of the differential framework through the introduction of the 
main definitions, results, and related examples. A detailed comparison with the existing literature 
is proposed in Section [Vj] Finally, LaSalle's invariance principle and horizontal contraction are 
presented in Sections IVIII and IVIIIL respectively. Conclusions follow. 
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II. Notation and preliminaries 

We present the differential framework on general manifolds by adopting the notation used in 
[DO and (HI. Given a manifold Ai of dimension d and a point x of Ai, we denote the tangent 
space of Ai at x by T x Ai. The tangent bundle of Ai is denoted by TAf. 

Given two smooth manifolds Ai\ and Ai 2 of dimension d\ and d 2 respectively, consider a 
function F : A4i — > Ai 2 and a point x E All, and consider two charts </?i : 14 C A(i — » R dl 
and 9? 2 : ZYp^) C A^2 - > R d2 defined on neighborhoods of x and We say that F is of 

class C k , k G N, if the function F = (p 2 o F o ip^ 1 : M dl — >■ ]R d 2 is of class C fc . We say that F 
is smooth (i.e. of class C 00 ) if F is smooth. The differential of F at x is denoted by DF(x)[] : 
T x A^i-)-r F ( a; )A / i2- h maps each tangent vector Sx E T x Mi to DF(x)[5x] E T F ^AiM. 

A curve 7 on a given manifold Ai, is a mapping 7 : / C R — >■ Al. A regular curve satisfies 
Z>y(s)[l] 7^ for each s E I. For simplicity we sometime use 7(3) or to denote Z>y(s)[l]. 
Following IfTOl Appendix A], given a C 1 and time varying vector field f : R x — > TA^ which 
assigns to each point x E M. a tangent vector /(t, x) G T^Al at time t, a C 1 curve 7 : / — > Ai 
is an integral curve of / if L>7(t)[l] = f(t, j(t)) for each t E L We say that a curve 7 : / — y Ai 
is a solution to the differential equation 7 = f(t, 7) on if 7 is an integral curve of /. 

Throughout the paper we will adopt the following notation. I n denotes the identity matrix of 
dimension n. Given a vector v, v T denotes the transpose vector of v. The span of a set of vectors 
{v u . . . , v n } is given by Span({«i, . . . , v n }) := {v \ 3A X , . . . A n G R s.t. v = YTi=i ^ v i}- Given 
a constant c G R we write R> c to denote the subset of [c, 00) C R. A locally Lipschitz function 
a : R> — y R>o is said to belong to class JC if it is strictly increasing and ct(0) = 0; it belongs 
to class /Coo if, moreover, lim r ^ +00 a(r) = +00. A function j3 : R>o x R>o — y R>o is said to 
belong to class ICC if (i) for each t > 0, /?(-, £) is a K function, and (ii) for each s > 0, f3(s, •) is 
nonincreasing and lim^oo f3(s, t) = 0. A distance (or metric) d : Al x Al — y R>o on a manifold 
Ai is a positive function that satisfies c?(x, ?/) = if and only if x — y, for each x, y G Ai and 
<i(x, 2) < d(x, y) + z) for each x, y, z E Ai. Given a set S C Af we say that S is bounded 
if swp x yeS d(x, y) < 00 for any given distance d on Ai. The distance between a set 5 and a 

1 We underline that the syntax DF(x)[v] follows the intuitive meaning of Differential of a function F : Mi — > A^2, computed 
at x £ and applied to the tangent vector v £ T^AI. D-F(a;)[t)] is replaced by more compact expressions like dF x v or 
in many textbooks. However, we found that the adopted notation makes the calculations more readable because of the clear 
distinction of the three elements F, x and v. 
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point x is given by d(A, x) := sup yeA d(y,x). Given a manifold M, a set of isolated points 
Vt C M. satisfies the following property: for any distance function d on Ai and any given pair 
xi,X2 in n, there exists an e > such that d(xi,x 2 ) > e. Given two functions / : Z — > y and 
g : X — >• Z, the composition fog assigns to each each p G <Y the value fog(p) = f(g(p)) G 3^- 
Given a function / : IR n — > W" 1 where we denote the (matrix of) partial derivatives by and 
we write ^p-^ x _ y for the partial derivatives computed at y G W 1 . 

III. Incremental stability and contraction 
Given the differential equation 

x = f(t,x) (1) 

where / is C 1 mapping, we denote by ip to (-,Xo) the solution to ([T]) from the initial condition 
x G M at time t » that is, ^0(^0,^0) — ^o- Following [33], we simplify the exposition by 
considering forward invariant and connected subsets C C M. for (OQ) such that ?p to (-,xo) is 
forward complete for every x G C, that is, ipt (t, xq) G C for each t > t - We assume that every 
two points in C can be connected by a smooth curve. 

The following definition characterizes several notions of incremental stability: 
Definition 1: Consider the differential equation (OQ) on a given manifold M.. Let C C Ai be 
a forward invariant set. The system (OQ) is 

(ZS) incrementally stable on C if there exists a distance d:MxjM-^ M> and a /C function 
a such that for any x\, x 2 G C and any £ > t > 

rf(^t (t,xi),^ (t,x 2 )) <a(d(x 1 ,x 2 )); (2) 

(7AS) incrementally asymptotically stable on C if it is incrementally stable and for any xi, x 2 G C 

lim d(if)t (t,xi),if) to (t,x 2 )) = 0; (3) 

t— >oo 

(ffiSJ incrementally exponentially stable on C if there exist a distance d, K > 1, and A > 
such that for any x\, x 2 G C and any t > t , 

d(ijH {t,x 1 ),ilH {t,X2)) < Ke' x ^ to) d(x 1 ,x 2 ). (4) 

j 

These definitions are incremental versions of classical notions of stability, asymptotic stability 
and exponential stability ITT31 Definition 4.4], and they reduce to those notions when either x\ or 



August 15, 2012 



DRAFT 



6 



x 2 is an equilibrium of (CO)- Global, regional, and local notions of of stability are specified through 
the definition of the set C. For example, we say that (QQ) is incrementally globally asymptotically 
stable when C = Ai. 

For At =M. n and for distances induced by norms on W 1 , the notions of incremental stability 
and incremental asymptotic stability given above are equivalent to the notions of incremental 
stability and attractive incremental stability of [16, Definition 6.22], respectively. For C = M n , 
the notion of incremental asymptotic stability is weaker than the notion of incremental global 
asymptotic stability of (31 Definition 2.1], since the latter requires uniform attractivity. 

Incremental stability of a dynamical system has been previously characterized by a suitable 
extension of Lyapunov theory fl3j . For Ai = M. n , the existence of a Lyapunov function decreasing 
along any pair of solutions is a sufficient condition for incremental stability [16, Theorem 6.30]. 
The key fact is in recognizing the equivalence between the incremental stability of x = f(t,x), 
x E W 1 , and the stability of the set A := {(x 1 ,x 2 ) E R 2n | X\ = x 2 } for the extended system 
X\ = f(t, X\), x 2 = f(t, x 2 ). As a direct consequence, incremental asymptotic stability is inferred 
from the existence of a Lyapunov function V(xi,x 2 ) for the set A with (uniformly) negative 
derivative along the dynamics f(t,xi), f(t,x 2 ), for any pair x\,x 2 . The extension to general 
manifolds is immediate. 

IV. Finsler-Lyapunov functions 

This section introduces a concept of Lyapunov function in the tangent bundle TAi of a 
manifold Ai. 

Definition 2: Consider a manifold Ai and a set of isolated points Q C M. A function V : 

M x TM — y K>o that maps every (x, 5x) E M. x T X J\A to V{x, 5x) E IR>o is a candidate 
Finsler-Lyapunov function for (0Q) if 

(i) V is of class C 1 for each x E Ai and 5x E T X M \ {0}; 

(ii) V(x, 5x) > for each x E M \ Q and each 5x E T X M \ {0}. 

(iii) for some p > 1, V(x, X5x) = X p V(x,5x) for each A > 0, x E Ai, and each 5x E T x Ai 
(homogeneity); 

(iv) for some p > 1, V(x, 8x\ + 5x 2 )p < V(x,5xi)p + V(x,5x 2 )p for each x E Ai \ f2 and 
each 5xi, 5x 2 E T x Ai such that 8x\ ^ XSx 2 and Ael (strict convexity). j 

For each x E Ai \fi, V is a measure of the length of each tangent vector 5x E T x Ai. The reason 
to call such a function V a "Finsler-Lyapunov function" is that it combines the properties of a 
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Lyapunov function and of a Finsler structure. The connection with classical Lyapunov functions is 
at methodological level: a candidate Finsler-Lyapunov function V is an abstraction on the system 
tangent bundle TAi, used to characterize the asymptotic behavior of the system trajectories 
by looking directly at the vector field f(t,x). Items (i) (regularity) and (ii) (positiveness) in 
Definition [2] are classical properties of a Lyapunov function. The regularity assumption (i) can 
be relaxed as in classical Lyapunov theory, see Remark [2] below. 

The connection with Finsler structures is immediate: for Q = and p — 1, V is a Finsler 



structure on Ai 11351 . while, for p > 1, the normalized function is a Finsler structure. There 
is no loss of generality in considering Finsler-Lyapunov functions which are homogeneous of 
degree p = 1 only. However, p greater than one may simplify the calculations in the application 
of the theory. Items (iii) and (iv) in Definition [2] are needed to turn the metric into a distance. 
In particular, through normalization, (iii) guarantees that the length of any curve 7 induced by 
V is independent on orientation-preserving reparameterizations of 7, while the combination of 
(iii) and (iv) guarantees that V(x, •) is a (possibly asymmetric) Minkowski norm in each tangent 
space. 

The Finsler structure of candidate Finsler-Lyapunov function V is a key property for the de- 
duction of incremental stability. In the next section, we will show that the incremental stability of 
the system can be inferred from the analysis of an infinitesimal approximation of the trajectories 
of the system - a local analysis - by the fact that V induces a well-defined distance on Ai, via 
integration. Following [[5l p. 145], 

Definition 3: [Finsler distance] Consider a candidate Finsler-Lyapunov function V on the 
manifold Ai. For any subset C C Ai and any two points 27, £2 £ Ai, let T(xi,X2) be the 
collection of piecewise C 1 , equally oriented curves 7 : / — > C, I := {s G ffi. | < s < 1}, 
7(0) = xi, and 7(1) = x 2 . 

The distance (or metric) d:MxjV(-> M> induced by V is given by 



r(x 1 ,x 2 )Jj 

where p is the degree of homogeneity of V(x, v) w.r.t v (see item (iii) in Definition [2]). j 
We consider curves whose domain is restricted to < s < 1 because the distance induced 
by V is independent from any orientation-preserving reparameterization of curves. With abuse 
of notation, in © we write 7(3) = Z>y(s)[l] to denote the directional derivative of a given 
piecewise C 1 function 7 at s, implicitly assuming that the differential is computed only where 




(5) 
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the function is differentiable. Points of non-differentiability characterize a set of measure zero, 
which can be neglected at integration. It is for the same reason that Definition [2] allows for 
functions V possibly equal to zero on isolated points in f2 C M.: this property simplifies the 
application of the main theorem in practical situations (See e.g. Example [3]), without affecting 
the integration. 

Example 1: We review specific classes of candidate Finsler-Lyapunov functions and classical 
distance functions. Consider C = M. = W l (for simplicity) and consider the Riemannian structure 
(5x 1 ,5x 2 ) x '■= 5xJ P(x)5x 2 for each x e M and each 8xi,Sx 2 G T x Ai, where P(x) is a 
symmetric and positive definite matrix in M nxn for each x G M. Then, the function V : A4 x 
TM — >■ IR>o given by V(x,5x) := (5x,5x) x satisfies the conditions of Definition [2] Moreover, 
from Definition [3l the distance induced by V is given by the length of the geodesic connecting 
x\ and x 2 . 

For the particular selection P(x) = I, V(x,8x) reduces to \5x\l- Thus, by normalization, 
d{x\,x 2 ) = Jq ds where 7 is the straight line 7(3) := (1 — s)xi + sx 2 . Therefore, 

d(xi,x 2 ) = Jq \x2 — xi\ 2 ds = \xi — x 2 \ 2 . Note that to induce distances d given by fc-norms 
d(xi,x 2 ) ■= \xi — x 2 \k, where k G N, k ^ 2, and Xi,x 2 G A4, a quadratic Finsler-Lyapunov 
function (i.e. a Riemannian structure) is too restrictive. Nevertheless, taking V(x, Sx) := \5x\k, 
we have that d(x\, x 2 ) = \x\ — x 2 \j-. j 

Example 2: We illustrate the necessity of conditions (iii) and (iv) for a candidate Finsler- 
Lyapunov function in Definition [2] for © to induce a distance. As a first example, consider the 
manifold M. = M 2 and take V(x, 5x) = 1 for each 1GK 2 and 5x G M 2 . V satisfies (i) and (ii) of 
Definition |2l while both (iii) and (iv) are not satisfied. Given any to points x%, x 2 G Ai, consider 
a generic curve 7 : I C M> — > M., I = [0, 1], such that 7(0) = Xi and 7(1) = x 2 . For any given 
normalizing exponent p, f I V( r y(s),'j(s))p = f T Ids = 1. Consider now a reparameterization of 
7 given by j k : I k -»■ M, h = [0, J], such that 7 fe (0) = x\ and 7 fc (|) = x 2 for any fe > 1. 
From Definition [3l the induced function d satisfies d(xi,x 2 ) < lim f T Ids = lim | = 0, 
for any given xi,X2 G .M. Thus, d is non-negative and satisfies the triangle inequality but 
d(xi,x 2 ) = for x\ 7^ x 2 . Therefore, d is not a distance. Note that a similar argument extends 
to V(x, Sx) = W(x) where W(x) is a positive and continuously differentiable function. 

As a second example, consider the simplified setting M. = M. Given the points and 
1, consider the curve 7fc(s) : [0, -] — > R such that 7jt(s) = ks, k G N>i. The function 
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V(x, 8x) := \5x\ Pl + \5x\ P2 is a candidate Finsler-Lyapunov function only if p x = p 2 , otherwise 
the homogeneity property in (iii) of Definition [2] does not hold. For non-homogeneous V, the 
integral in the right-hand side of © depends on the particular parameterization of the curve, for 
any selection of the normalizing exponent p G N. In fact, for any given p, and any given p 1 and 
p 2 , we have that V ( r yk( s ) , ik( s ))^ — Jq (k Pl + k P2 )p ds = | (k Pl + k P2 )p which preserves a 
constant value for any given k only when p = p\ = p 2 . j 
The reader will notice that the distance d induced by a candidate Finsler-Lyapunov function is 
not symmetric in general, that is, we may have d(x, y) ^ d(y, x) for some x, y G Ai. To induce a 
symmetric distance, it is sufficient to strength (iii) in Definition [2] to (iii)?, V(x, \v) = \X\ p V(x, v ) 
for each A, p > 1, x G Ai, and each v G T x Ai (absolute homogeneity, [35]). Note that 
adopting (iii) b reduces the generality of the class of Finsler-Lyapunov functions excluding, for 
example, Randers metrics (U Section 1.3]. Finally, when coordinates are available, conditions 
(iv) of Definition [2] can be replaced by the (possibly simpler) condition that, for each for each 

2 

x G Ai \ fl, the Hessian of V(x, -)p (p is the degree of homogeneity with respect to 5x) is 
positive definite at every point v of T x Ai \ {0}, 11351 Condition F.iii]. 

V. A Finsler-Lyapunov theorem for contraction analysis 

For any given candidate Finsler-Lyapunov function V : xTAI — > M>o, we use DV(x, 5x) [■} 
to denote the differential of V computed at (x, 5x) 6Mx T x Ai. Note that Ai x TAi itself has 
the structure of a manifold. Thus, DV(x, Sx)[-] is a mapping from the tangent space TL t $ x )(Ai x 
TAi), whose tangent vectors are represented as pairs (v 1 , v 2 ), to the tangent space T V (x,Sx)^>o — 
E, that is, (vi,v 2 ) t-> DV(x, 5x)[v 1 ,v 2 \. 

Theorem 1: Consider the system (OQ) on a smooth manifold Ai with / of class C 1 , a connected 
and forward invariant set C, a set of isolated points OCM, and a function a : M> — > IR>o- 
Let V be a candidate Finsler-Lyapunov function such that 

DV(x,8x)[f(t,x), Df(t, x)[0,5x\] < -a(V(x,6x)) (6) 

for each t > t , each x G C C Ai, and each 5x G T x Ai. Then, dU) is 
incrementally stable on C if a(s) = for each s > 0; 
(/AS) incrementally asymptotically stable on C if a is a /C function; 
(.ffiS) incrementally exponentially stable on C if a(s) = Xs > for each s > 0. 
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The conditions of the theorem for incremental stability are reminiscent of classical Lyapunov 
conditions for stability, asymptotic stability and exponential stability lfl"5l Chapter 4], extended 
to the tangent bundle TAi. We say that the system (OQ) contracts V in C if © is satisfied for 
some function a of class /C. V is called the contraction measure, and C the contraction region. 

From the proof of the theorem it will be clear that © enforces the decay of the function 
V along each solution vfj to (t,x ), x G C. Thus, exploiting the relation between V and Finsler 
structure, the contraction of the structure along if} to (t,x ) (locally - in each tangent space) guar- 
antees, via integration, that the distance between any pair of solutions ip to (t,xi) and ^t (t,x 2 ), 
Xi, x 2 G C, shrinks to zero as t goes to infinity. A graphical illustration is provided in Figured] 



Figure 1. A graphical illustration of the contraction of the distance induced by Condition {6]l on the solutions to (Q3. The 
Finsler-Lyapunov function assigns a positive value to each pair (7(5), 7(s)). The length of the curve 7 is given by the integral 
of the Finsler-Lyapunov function along 7, represented by the shaded area. 



The incremental Lyapunov approach proposed in [3], establishes incremental stability by 
checking a pointwise geometric condition in the product space Ai x Ai . In contrast, the differen- 
tial approach proposed here establishes incremental stability by checking a pointwise geometric 
condition in the tangent boundle TAi. Several earlier works have adopted this approach in a 
Riemannian framework, focusing on quadratic functions V(x,5x) = 8x T P(x)Sx in Euclidean 
spaces ( lf22l . Slotine's contraction lfT8l Section 3.4], [|2l Appendix II], [28, Equations (13), 
(14), and Lemma 2.2]). There are a number of reasons to consider Finsler generalizations of 
Riemannian structures for contraction analysis, some of which are illustrated in the next section, 



d(7(0), 7 (l)) 




d(7(0), 7 (s)) 
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where we report a detailed comparison between the conditions proposed in Theorem Q] and 
several results available in literature. 

Before entering into the details of the proof, we present a scalar example that illustrates the 
value of non constant Riemannian structures in nonlinear spaces. 

Example 3: For Ai = S 1 consider the dynamics 

# = /(#) : = -sin(tf). (7) 



2 



The tangent space at every point d G M. is given by R. The naive choice Vi(i?, 5t?) := 
corresponds to a flat Riemannian structure on 5 1 . Then, for any given compact set C C (— | , | ) 
© yields 

dVl ^' 5 ^ ( ^P-) W = - cos(^)^ 2 < -eVi(fc?) for each G C, (8) 
cw \ ov J 

where e > (sufficiently small). From Theorem [T] we conclude that © is incrementally 
exponentially stable on compact sets C C (— f,f) such that G C (to guarantee that C is 
forward invariant). For C = [— ~, ~] we have only incremental stability, since cos($) = at 
= f. From Definition [3l note that the distance induced by V\ is given by — t? 2 | - 
The true contracting region is captured with the choice V 2 '■ (S 1 \ {7r}) x R — > R> given by 
V 2 = tan (|) 2 5d 2 . Despite the identification of each TgAi with R, the measure of the "length" 
of 5§ given by V 2 now depends on particular the tangent space Note that V 2 satisfies each 

condition of Definition [2] with a set Vl = {0}. Moreover, V 2 is well defined in S 1 \ {n} since 
tan (|) 2 — > oo as — > ir. For any given compact set C C (S 1 \ {7r}) such that G C, using 
the identity tan(|) = ^ff^ , © yields 

^ m + ^(m S ^ = -2tan(f ) ^g T ^-2tan(f ) 2 cos(^ 



^ W J ~~~ " 2 / cos(i9)+l 

(9) 



-2 tan(f ) W - 2 tan(f ) 2 cos(tf)c^ 2 



= -2(1 + cos(tf))F 2 (tf, 
< -eF 2 (^, &?) for each G C, 

where £ > 0. Thus, by Theorem [fl © is incrementally exponentially stable on C. j 

Proof of Theorem\l} The proof is divided in four main steps. 
(i) Setup: Finsler structure and parameterized solution. 

For any two points Xi,x 2 G Ai, let T(xi,x 2 ) be the collection of piecewise C 1 , equally 
oriented curves 7 : / -> C C .M, / := {s G R | < s < 1}, connecting xi to x 2 , that is, 
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7(0) = X\ and 7(1) = x 2 . Considering the distance d induced by V (Definition [3]), we have that 

d(x u x 2 ) = inf fv{j{s),Dj{s)[l))pds. (10) 
r(xi,» 2 ) J j 

For any two initial conditions xi,x 2 E C and any given e > 0, consider now a regular C 1 
curve 7:/->CcM such that 7(0) = X\, 7(1) = x 2 , and ^ 

Jv(?y{s),I*y(s)[l])pds< {l + e)d{x 1 ,x 2 ). (11) 

Let ^ io (-,7(s)) be the solution to (QQ) from the initial condition 7(5), for s E I, at time i - 
Precisely, ^t (-,7(-)) ls a function from R x / to M.. Its differential at (t,s) is the mapping 
£tyto(M)H fr o m x /) = R 2 to T^ to(tMs)) M. A simple base for % s) (R x J) is given 

by the tangent vectors and ^ which are denoted respectively by (1,0) and (0,1). Thus, 
ift to (•,7(-)) satisfies ([U) in the sense that 

Dip to (t,j(s))[l,0) = f(t,i/j t0 (t,%s))) for each t > t and s E I. (12) 
Clearly ^t (to, 7(0) = 7(0 tnus ' fr° m (fTTj) . we have that 

^V(^ (t ,7(s)),^ io (to,7(s))[0,l])^s < (l + £ )rf(x!,x 2 ). (13) 

For each t > t and s E [0, 1] the differential of ip in the direction (1,0) characterizes the time 
derivative of the parameterized solution if} to (', 7 y(s)). Instead, the differential of ip in the direction 
(0, 1) E T( tjS )(R x [0, 1]) characterizes at each s the tangent vector to the curve ipt (t,^y(-)), for 
fixed time t. Following [TT8l . we call this tangent vector virtual displacement. Thus, combining 
integration of the displacement along (0, 1), time derivative along (1, 0), and ©, we can establish 
contraction of the distance in (PTOl) along the solutions to (Q~|). 

(ii) The displacement dynamics along the solution ^ (-,t"(s)). 

Consider the function Sijj to (-,-) : R x I — > TM. given by the tangent vector 8ip to {t,s) = 
D^ to (t, t"(s))[0, 1] for each t > t and s E I. We study the time derivative of 8^ to (t,s) by 

2 By using a generic curve 7 £ C 1 which satisfies dl It we do not need to assume the existence of geodesies and we simplify 
the exposition by avoiding the analysis of points of non-differentiability. 
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looking at the differential of Si[) to (-, •) at (t, s) in the direction (1, 0), that is, 

D6i> to (t,s)[l,0] = D(^ to (-,7(-))[0,l])(M)[l,0] (14a) 

= J D(^ to (-,7(-))[l ) 0])(t, S )[0,l] (14b) 

= Df(t,i> to (t,j(s)))[0,l] (14c) 

= Df(t,ij to (t,^(s)))[0,Di; to (tMs))[0,l}} (14d) 

= Df(t,^ to (t,l(s)))[0,8^ to (t,s)}. (14e) 



(114ab follows from the definition of 5ift to (-, •). (|14b|) follows from the fact that ip is a C 1 function. 
In coordinates we would have -^ipt {t, 7(s)) = ^i^o^ 7~( s )) for a C 1 function W CS> 
follows from (fl"2)) while (|14dl) follows from the chain rule. Because the differential in (|14cl) 
is a mapping from T (t)S )(lR x J) to T f{t ^ to{t ^ s))) TM (it is constructed from f(t, V> to (t, t(s))) 
as a function of t and s). In coordinates, (114cb reads -§^f(t, ^ (t, 7(s)))- The chain rule in 
(|14d| ) uses the combination of two differentials: T( 4)S )(R x /) — > Tu^ t ( t)! y( s )))(R x Ai) — > 

T f(t,^ (mm™- In coordinates & reads |x= ^_ (s)) ^ (*, ^(s))- Finall y> fl© fo 1 " 
lows from the definition of 6ip to (t, s). 

fiiij 77ze dynamics ofV along the solution i/; to (-,j(s)). 

Consider the function V:lx/-^ R> given by V(t, s) = V(ijj to (t, j(s)), 5ip to (t, s)) for each 
t > t and s e /. Note that V has a well-defined time derivative -^V(t, s) since V(t, s) G R>o 
for each t and s. 

jV{t,s) = DV(t,s)[l,0} (15a) 

= W(^,7(s)),^M))[^ (15b) 

= ^(.^[/(^(^s))),])/^^^)))^^^)]] (15c) 

< -a(K(^ (t,7(s)),# 4o (t,s))) (15d) 

= -a(F(*,s)). (15e) 

(|15b| ) follows from the application of the chain rule, while (|15c|) follows from (PT21) and (PT4l) . 
( |15d| ) is enforced by ©. Finally, (|15eb follows from the definition of V Note that (|15bl) in 

" <9^ (t,7(s)) " 



coordinates is given by the product 



dV(x,8x) dV(x,Sx) 
dx dSx 



\(x,8x)={ipt (i,7(s)),5^t (t,s)) 



The coordinate representation for (I15cl) is similar. 



a* 
Ft 
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(iv) Incremental stability properties. 

(IS) Incremental stability: if a(s) — for each s > then 

V(t, s) < V(t , s) for all t > t and s e I. (16) 

Therefore, for each t >t , 

d(ipto{t,x 1 ),ip to (t,x 2 )) < JV(t,s)pds = Jv(t ,s)vds < (1 + e)d(x h x 2 ) (17) 

i 

where the first inequality follows from the definition of induced distance in ([5]), since V(t, •)» 
is constructed on a "particular" curve connecting ?p to (t,xi) and tp to (t,X2); the identity follows 
from (fT6l ); and the last inequality follows from (fT3l) . 

(IAS) Incremental asymptotic stability: if a is a /C function then (IS) holds, and by the 
comparison lemma [121 Lemma 3.4] and Q~51 Lemma 4.4], there exists a ICC function (5 such 
that 

V(t, s) < /3(V(to, s),t- t ) for all t > t Q and s E I. (18) 
Therefore, for each t >t , 

d{ip t0 {t,x 1 ),ip tQ {t,x 2 )) < Jv{t,s)Us < J /3{V{to,s),t-t fpds (19) 
from which we get 

limtZ(V^(t,a;i),Vto(*^2)) < lim / (3(V(t ,s),t -t )*ds = 0. (20) 

t— too t— >oo J j 

(IES) Incremental exponential stability: if as = As > for each s > then, by [fT5l Lemma 
3.4], we get 

V(t, s) < e- x{t - to) V(t , s) for all t > t and s e /. (21) 
Therefore, for each t > t , 

d(ij t0 (t,x 1 ),ij t0 (t,x 2 )) < jv{t,s)Us < e-^ to) Jv(t ,s)Us < (l+e)e-^ t - to) d(x 1 ,x 2 ). 

(22) 
■ 

The proof of Theorem \T\ generalizes the argument proposed in the proof of [|33l Lemma 1] 
and [27, Theorem 5] to general manifolds and Finsler structures (the proof provided in ll2~7ll is 
developed for Euclidean spaces using matrix measures). An equivalent proof to Theorem Q] for 
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incremental exponential stability and V restricted to Riemannian structures can be found in |[2l 
Appendix II]. 

Remark 1: From (fTTT) and (fTTT) . for any given converging sequence Ek G M>o> li m £fc = 0, 

k— >co 

we can construct a sequence of curves 7^ : I k — > Ai such that 

lim / V(j k (s),Dj k (s)[l))pds < lim (1 +e k )d(x 1 ,x 2 ) = d(x 1 ,x 2 ). (23) 

fc— >oo fc— >oo 

Therefore, in the limit of A; — )■ 00, the conditions of Theorem \T\ guarantees incremental stability 
with the stronger property that 

d(ijjt (t,x 1 ) } ip to (t J x 2 )) < d(xi,x 2 ) Vt>t ,Vxi,x 2 eM, (24) 

which becomes a strict inequality for t > t when © is satisfied by a K, function a, as shown 
by ©. j 
Remark 2: The result of Theorem Q] can be extended to piecewise continuously differentiable 
and locally Lipschitz candidate Finsler-Lyapunov functions V (relaxation of Item (i) in Definition 
O. The key observation is that the decrease of the distance between any two solutions is preserved 
also if (TT3T) holds for almost every t and s. With this aim, let V C Ai x TAi be the set of 
nondifferentiability points of V. Under the hypothesis that © holds on every point where V 
is differentiable, if for each t > t , and each (x, Sx) E V such that Sx G T x Ai, the pair 
(f(t, x), Df(t, x)[0, Sx]) does not belong to the tangent cone to V at (x,5x), then (fT5l) holds 
for almost every t and s . In fact, this geometrical condition guarantees that for any given 
solution ijj to such that (ip to (t, x), Dtjj to (t, x)[0, 5x]) G V, there exists e > which guarantees 

(Vfo (t, ar) , .D Vi (t» [°> ^ v for ever y r G (t, t + e] . j 
We conclude the section by emphasizing the geometric (or coordinate-free) nature of Theorem 
\T\ Specifically, consider a candidate Finsler-Lyapunov function V x expressed in a given set of 
coordinates X C M. d Ai, and suppose that V x satisfies the conditions of Theorem Q] for 
the system x = f(x). Consider also a coordinate transformation h : W 1 — > W 1 that smoothly 
maps every point x G X to y = h(x) G y C M. d . The dynamics of y is thus given by 
V — g{y) := Dh(h~ l (y))[f(h~ l (y))}. Then, in the new set of coordinates, the function V y given 
by V y (y,8y) := V x (h^ 1 (y), Dh~ 1 (y)[5y}) satisfies the condition of Theorem Q] for y = g(y). 
This follows from the fact that both V y and V x are coordinate representations of the same 
"coordinate-free" object V : M. x TAi — >■ M>o which assigns to each point (p, <5p) the value 
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V(p,5p) := V x (^ x 1 (p),Dcp x 1 (p)[5p}) or, equivalently, V(p,5p) := ^(^/(p), D<p y x (p)[5p]) 
where tp y := (p x o h~ l and ^ M. 

VI. Revisiting some literature on contraction 
A. Riemannian contraction, matrix measure contraction, and incremental stability 

For a historical perspective on contraction the reader is referred to [fl"3"1 . and related concepts 
in [|22| and [|33l . We propose here a detailed comparison with selected references from the 
literature. First, we consider results on contraction based on matrix measures ll27l . Il33l and matrix 
inequalities [23 J. We recast these results within the differential framework proposed in Theorem 
[Q by suitable definitions of state-independent Finsler-Lyapunov functions V(x,5x). Then, we 
consider results based on Riemannian structures [fT8l . 0, and we show that they coincide with 
the (IES) condition of Theorem \T\ for a function V(x,5x) defined by the Riemannian structure. 

The reader will notice that these two groups of results are essentially disjoint. The equivalence 
between the conditions based on matrix measures and the conditions based on Riemannian 



structures can be established only for quadratic vector norms \x\p — Vx T Px or, equivalently, 
for state-independent Riemannian structures (Sx, 5x) = 5x T PSx. However, both groups of results 
fall within the proposed differential Finsler-Lyapunov framework. 

Since the approaches in [1271 , 11331 , [|23l consider an Euclidean space, as a first step we rewrite 
© for Ai = M. n . Consider the system x = f(t,x), x E W 1 , and denote the Jacobian of / 
at x by J(t,x) := Df(t,x)[0,-] = 9 ^q^ ■ Then, for a candidate Finsler-Lyapunov function 
V : R n x W n -> E> , © can be written as 

^lM/( t , x) + jfc x)6x < - a(V ( x , Sx) ). (25) 

OX oox 

The approach proposed in ll27ll and [|33l is based on the matrix measure of the Jacobian 



J(t,x). For instance, given a vector norm | ■ | in IR n and its induced matrix norm, the induced 
matrix measure \l of a matrix A G R nxn is given by p(A) := lim — — -, [37, Section 



3.2]. Then, following 11331 Definition 1 and Theorem 1], let C be a convex set, forward invariant 
for the system x = f(t,x). f is a C 1 function. If 

fi(J(t, x)) < — c < for each x G C and each t > 0, (26) 

then the system is incrementally exponentially stable with a distance given by d(xi,X2) = 
\x\ — £2 1- Moreover, by ll33l Lemma 4], the same result hold for non convex sets C that satisfy 
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a mild regularity assumption, and it guarantees incremental exponential stability with a distance 

function d(xi, 22) < K\x\ — x%\ for some K > 1. 

Condition (|26l) is equivalent to (1251) for the state-independent Finsler-Lyapunov function V 

given by V(x, 5x) = \5x\ and a(s) = cs. This follows from 

dV(x, 8x) V(5x + hJ(t,x)5x) — V(x,5x) 

— J(t,x)ox = hm 

00 x h^0+ h 

\I + hJ(t,x)\\5x\ - \8x\ 



< lim 



h-+0+ h 

= lim l J + -V(,,fa) < 27 > 

h->-0+ h 

= fi(J(t,x))V(x,8x) 

= -cV(x, Sx) for each t > 0, x G C, 5x e W 1 . 
The approach proposed in ll23l (and in [|39l Chapter 5, Section 5] for time-invariant systems) 
use matrix inequalities based on the Jacobian J(t, x) and on two positive definite and symmetric 
matrices P and Q. The result ll23l Theorem 1] is equivalent to incremental exponential stability of 
Theorem Q] for V is restricted to the constant Riemannian structure 5x T P5x. Consider the system 
x = f(x,w(t)) where / is a C 1 function and w : M>o — > W C M m is a C 1 exogenous signal. 
Thus, f(x,w(t)) is a time-varying C 1 function. Applying Theorem Q] to V(x,5x) = 5x T P5x, 
incremental exponential stability holds if 

dV(x 5x)df(x,w) Sx = ^ T f p df(x,w) + df(x,w) T p \ Sx < _ xv{x ^ 5x) = „ XSx t PSx 
oox ox \ ox ox I 

(28) 

for some A > and for every Sx E W 1 and w E W. The right-hand side of (1281) can be replaced 
by —5x T Q5x, for some matrix Q = Q T > (for any given Q, we can always find A sufficiently 
small to guarantee Q > \P, and vice versa). Therefore, the condition in (1281) is equivalent to 
the existence of positive definite and symmetric matrices P and Q such that 

P df ^ w) + d J^l T p < -Q (29) 

OX OX 

which is 11231 Eq. (8), Theorem 1]. The induced distance is the quadratic form d(xi,X2) = 



a/ (xi— x 2 ) T P(xi — x 2 ). See also [|24| and Section IVI-BI in the present paper. 

Conditions for contraction based on quadratic structures 5x T M(x)5x are provided in the 
seminal contraction paper lfT8l (we consider the time-invariant case only). lfT8l Definition 2 and 
Theorem 2] establish incremental exponential stability for x = f(t,x) by requiring, using the 
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notation of [ 183], that the inequality 

6x T (j(t, x) T M(x) + M{x)J{t, x) + M(x)) 5x < -X5x T M(x)5x (30) 

is satisfied for every x and 5x, for some A > 0. Note that Sx T M(x)5x is a short notation 
for jj^(5x T M{x)8x)f{x). Therefore, taking V(x,5x) = 5x T M{x)8x, the relation between ( 1301 ) 
and (1251 ) for incremental exponential stability is immediate. The same argument illustrates the 
relation between the differential approach proposed here and the results in [2, Appendix II] and 
HOI Definition 2.4 and Theorem 2.5] (for this last paper, the differential equation x = f(x,u), 
where u is an input signal, is casted to the form (OQ) by considering the time-varying vector field 
J(t,x) := f(x,u(t)))- 

We conclude the section by considering the incremental Lyapunov approach in Q, [1261 . The 
key observation is given by (31 Lemma 2.3 and Remark 2.4] and [26, Appendix A.l] which shows 
the equivalence between the incremental stability of x = f(t,x), x G M n , and the stability of 
the set A := {(xi,x 2 ) G M 2n | X\ = x 2 } for the extended system x\ = f(t,Xi), x 2 = f(t,x 2 ). 
Thus, to show asymptotic stability of the set A, a Lyapunov function V(xi, x 2 ) must be positive 
everywhere but on A, that is 

a(\xi - x 2 \) < V(xi,x 2 ) < a(\xi - x 2 \), (31) 

for some a, a G /C; and the derivative of V(xi,x 2 ) along the solutions of the system must 
decrease for xi,x 2 A, which is established by enforcing 

dV(xi,x 2 ) . , dV{x 1 ,x 2 ) \ ~ / j h 

— f{t, xi) H — f(t, x 2 ) < -a(\xi -x 2 \) (32) 

for each pair Xi,x 2 G M n , where a G K. Indeed, an incremental Lyapunov function is essentially 
a Lyapunov function for the extended system which measures directly the distance between any 
two points x\ and x 2 . 

The differential framework proposed here does not use a Lyapunov function to study directly 
the time evolution of the distance between any two solutions. Instead, a lifted Lyapunov function 
on the tangent bundle is used to characterize the contraction of the infinitesimal neighborhood 
of each point x - a local property - to infer indirectly the contraction of the distance - a global 
property - via integration. Applications suggest that it can be considerably more difficult to 
construct a distance than the associated differential structure. 
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B. Contractive systems forget initial conditions 

All the solutions of a contractive system converge to a unique steady-state solution. This 
feature is exploited in control design ll38l . Il23l . Il24ll . Ifl4l . for example in tracking, by inducing 
an attractive desired steady-state solution via the feedforward action of exogenous signals (that 
preserve the contraction property), or in observer design, by a suitable injection of the measured 
output. In what follows we revisit these results, showing that a particular application of Theorem 
Q] entails the sufficient conditions for convergent systems in 11231 , Il24ll . and we formulate a 
proposition whose conditions parallels the relaxed contraction analysis proposed by 113811 , lfT4ll . 
through the notion of virtual system. 

Following [|23l and Il24ll . consider the system x = f(x, w(t)) where w is an exogenous signal. 
Define f(t,x) := f(x,w(t)) and suppose that Theorem Q] holds for x = f(t,x). Then, by 
incremental asymptotic stability, the solutions of the system converge towards each other, thus 
every solution converges to a steady state solution x*(t) = f(x*(t),w(t)) induced by w. This 
results parallels [23, Property 3]. In particular, Theorem Q] applied to x = f(t,x) = f(x,w(t)) 
recovers ll23l Property 3] when V = Sz T PSz (flat metric) and a(s) = —ks, k > 0. 

Following ll38l and |[T4ll . consider the system (OQ) given by x = f(t,x) and a new system of 
equations 

z = f(t,z,x) such that f(t,x,x) = f(t,x),feC 1 . (33) 

(1331 is the so-called virtual system, ||38l . (1331 arises naturally in tracking and state estimation 
problems where, possibly, CD is the reference system and the controlled/observer system is given 
by (|33l) . For example, f(t, z, x) = f(t, z) + K{z — x) may represent a tracking controlled system 
with state-feedback K[z — x), while f(t, z, x) = f(t, z) + L(y z — y x ) may represent an observer 
dynamics with output injection L(y z — y x ). Inspired by 11381 and |[T4l . we provide the following 
proposition, a straightforward application of Theorem \T\ 

Proposition 1: Consider the system (OQ) on a smooth manifold Ai with / of class C 1 , and a 
connected and forward invariant set C x C M. for (OQ). Consider ([33]) and suppose that the set 
C z C M. is connected and forward invariant for (1331) . Given a K, function a, let V be a candidate 
Finsler-Lyapunov function for d33l (Definition [2]) such that 

DV(z,Sz)[f(t,z,x),Df{t,z,x)[0,Sz,0]] < -a(V(z,Sz)) (34) 

for each t > t , each x E C x (uniformly in x), each z E C z C Ai, and each Sz E T z Ai. Then, 
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for any given initial condition x G C x , and any initial condition z G C z , each solution ip^(t, z ) 
to (1331) converges asymptotically to the solution cpf (t, Xo) to (QQ). j 
Combining the virtual system decomposition (|33T) with Proposition \T\ is useful for applications 
like tracking and state estimation, but also as an analysis tool. In fact, if Proposition Q] holds and 
(1331) converges to a given steady-state solution z* uniformly in x, then all solutions of (OQ) converge 
to that solution. The conclusion of Proposition \T\ is a consequence of Theorem [Q considering 
the solution (p^ Q (t,Xo) to (0Q) from a given initial condition x £ £e, the dynamics (|33l) can be 
rewritten as the time-varying dynamics i = f(t,z) := f(t,z,(pt (t,x )), and (l34l) guarantees 
that the conditions for incremental asymptotic stability of Theorem \T\ applied to z — f(t, z) are 
satisfied. Therefore, for any given initial conditions Zi,z 2 , the solutions ip^ Q (t, zi) and ip^ (t,z 2 ) 
converges towards each other, that is, lim d(ipf (t, z{), ift? (t, z 2 )) = 0. The conclusion of the 

t— s-oo 

proposition follows by noticing that when z 2 = x , we have that ^ Q (t,z 2 ) = 4>t {tj x o) (since 
f(t, x, x) = f(t, x)). Thus, from every initial condition z\ G C z , lim d{ipi ' (t, zi), ipf (t, x )) = 0. 

t— too 

Similar conditions are provided in [|38l and IfPfll for Pdemannian metrics V(z, Sz) = 5z T P(z)Sz. 

VII. LaSalle-like relaxations 

A very first step of Lyapunov theory is to relax the strict decay of Lyapunov functions by 
exploiting the invariance of limit sets. We show that this important relaxation readily extends 
to Finsler-Lyapunov functions. We develop the analysis for the particular case of time-invariant 
differential equations x = f(x). 

Theorem 2: [LaSalle invariance principle for contraction] Consider the system x = f(x) on 
a smooth manifold Ai with / of class C 1 , a continuous function a : M. x TAi — >■ M>o, a set 
of isolated points OcM, and a connected set C C Ai, forward invariant for x = f(x). Let V 
be a candidate Finsler-Lyapunov function such that 

DV(x,6x)[f(x),Df(x)[6x]] < -a(x,6x) (35) 
for each x G C C Ai, and each 5x G T x Ai. Define 

n T := {(x,5x) eMx (TM \ {0})\a(ip(r,x),Di/;(T,x)[0,8x]) = 0} (36) 

where ip(-,x) is the solution^ to x = f(x) from the initial condition x G Ai. Then, x = f(x) is 
incrementally asymptotically stable on C if there exists r > such that il fl II T = 0. j 

3 Without loss of generality we drop the subscript to- 
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Proof: Incremental stability follows from (1351 ) by Theorem [TJ For incremental asymptotic 
stability we have to show that lim f T V(ip(t, 7(s)), Dip(t, 7(s))[0, 1]) = 0, for any given curve 

t— >oo 

7 : I^-M. This is equivalent to establish that lim ^(-^(^(s)), Di/j(t,"f(s))[0,l\) = 0, for each 

t— >oo 

s e I. Suppose that there exists a T > such that V(ip(t, j{s)), Dip(t, 7(s))[0, 1]) > c> for 
each s£l and t >T. Without loss of generality we can consider T = Q (1351) and continuity 
of a imply that a(ip(t, 7(s)), Dip(t, 7(s))[0, 1]) = 0, for each £ > and s E I. Therefore, 
V^^, 7(s)), DiJ}(t, 7(s))[0, 1]) > c > for each s G J and i > implies the existence of a 
set G .M | a(ip(t, x), Dip(t,x)[0, Sx\) = for each t > 0}. But this set is empty since 

n n n T = for some r > 0. ■ 

To the best of author's knowledge, an invariance principle has not appeared in the literature 
on contraction. This illustrates the potential of a Lyapunov framework for contraction analysis. 

We illustrate the use of Theorem [2] in the following (linear) example, where we take advantage 
of classical observability conditions. Example 0] illustrates a general class of models in power 
electronics for which incremental tools are frequently used 11301 . 

Example 4: Consider the following averaged equations of a single-boost converter flU 

Lx L = 

Cx C = UX L - jjX C 

where xl is the inductor current, xc is the capacitor voltage, and E is the input voltage. The 
quantities L, C, and R are respectively the inductance, the capacitance and the (load) resistance 
of the circuit. 

We claim that for any given constant input u* ^ 0, and any constant positive value of the 
circuit quantities L, C and R, the system is incrementally asymptotically stable. Note that (1371 ) is 
a time-invariant linear system for u = it*, so that a natural candidate Finsler-Lyapunov function 
is provided by the incremental energy V(x, 5x) = ^(LSx^ + CSxq). In fact, 



-uxc + E 
i 



(37) 



DV(x,5x)[f(x,u*),Df(x,u*)[5x,0]} 



5x L 


T 


-u* 




5x L 


5x c 




U R 




5x c 



Sx 



c 



R 



<0, 



(38) 



where a(x,5x) = ^p. By d36l) . considering ip(t 7 x) = e At x, we have that I1 T := {(x, 8x) G 
M x TM | 5x T (e Ar ) T [g ?] e Ar 5x = 0}. Thus, for any given r > 0, we have that n H I1 T = 0. 



4 By considering the new solution ip(t, ip(T, j(s))), t > 0. 
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Incremental asymptotic stability follows from Theorem 2 (for the linear nature of the system, 
the incremental asymptotic stability is exponential). j 
Remark 3: For a time-varying differential equation ([I]), a possible formulation of invariance- 
like conditions for asymptotic stability is given by the inequality 

DV(x,6x)[f(t,x),Df(t,x)[0,6x\] < -a(t,x,8x)V(x,8x) (39) 

for each t > to, each x G C C Ai, and each 8x G T x Ai, where V is a candidate Finsler-Lyapunov 
and a : IR> xMx TM — > M> . Incremental asymptotic stability on C holds if 

lim / a(r, ■0* o ( r > x), Di/j to (r, x)[0, 5x\)dr = oo for each x G C and Sx G T X M. 

*- >0 ° An 



(40) 

In general, (1401) is established by relying on further analysis of the solutions of the system^. 

By Theorem \T\ (|39l guarantees incremental stability. To see why (|39l and (|40l ) guarantee 
incremental asymptotic stability, one has to follow the proof of Theorem \T\ up to Equation ( fT5T ), 
by replacing each quantity fe)) by a(t, x, 8x)V(x, 8x). From there, using the definition 

a(t,s) := ^0(^,7(5)), Z^^ to (t,7(s))[0, 1]), by comparison lemma [15, Lemma 3.4] we get 
V"(£, s) < e~ jt o alT ' SJ T V(t ,s) for all t > t and s 6 J, which combined with (1401) guarantees 
that 

lim d(ip to (t,Xx), i) to (t,x 2 )) < lim Je~p ^o a{T ' s)dT V(t ,s)pds 

t— too t— >oo 

< ( maxV(t ,s)v) lim f r e~p J ^ s)dT ds (41) 

= 0. 

VIII. Horizontal contraction 
A. Contraction and symmetries 

Theorem [H guarantees contraction among the solutions of a system in every possible direction. 
This result can be easily extended to capture contraction with respect to specific directions 
- a relevant feature for contraction analysis in presence of symmetries like, for example, in 
synchronization problems or observer design. 

The generalization of Theorem \T\ is based on the introduction of horizontal Finsler-Lyapunov 
functions on a manifold M, whose induced metric d is tailored to the particular problem of 

5 The differential Dipt (i+ r i x )[0, ■] assigns to each tangent vector Sx 6 Ty, t (t,x)M the tangent vector Dipt (t+T, x)[0, Sx] £ 
T^ t (t+r.x) M. Thus, Dtpt (t + t,x)[0,Sx] represents the evolution of the tangent vector Sx along the solution t/>t after r 
units of time. 
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interest. These functions are positive only on a suitably selected (horizontal) subspace % x C 
T x Ai, for each x G Ai, which characterize the set of directions (tangent vectors) taken into 
account by the Finsler structure. 

Definition 4: [Horizontal Finsler-Lyapunov function] Consider a manifold Ai of dimension d 
and a set of isolated points QcM. For each x G Ai, consider a subdivision of T x Ai into a 
vertical distribution V x C T x Ai given by 

V x = Span({t>i(a;), . . . ,v r (x)}), 0<r<d (42) 

where each Vi : Ai — > T x Ai, i = l,...,r, is a function of class C 1 (vector fields), and a 
horizontal distribution T-L x C T x Ai complementary to V x , i.e. V x © T-L x = T x Ai, given by 

% x = Span({/ii(x), . . . , h q (x)}), < q < d — r (43) 

where each hi : Ai — > T x Ai, i = 1, . . . , q, is a function of class C 1 . 

A function V : Ai x T.M -»■ IR> that maps every (x, for) G Ai x T X .M to V(x, fo) G M>o 
is a candidate horizontal Finsler-Lyapunov function for (OQ) on % x if 

(i) V is of class C 1 for each x £ Ai and 5x £H X \ {0}; 

(ii) V^(x, <5x) > for each x £ Ai\£l and 5x G "% x \ {0}. Moreover, V(x, <5x) = V^x, 5xh) for 
each igM and 5x G T^-A^ \ {0} where 5a; = <5x/j + 5x v , 5xh G and fo^ G V x ; 

(iii) for some p > 1, V(a;, A5x) = A p V(x, 5x) for each A > 0, x G .M, and <5x G 

ill 

(iv) for some p > 1, F(x, <5xi + 5x 2 ) p < l/(x,&ri)p + K(x,5x 2 ) p for each x G .M \ f2 and 

Sxi, 6x2 G \ {0} such that 5x\ ^ A#x 2 and A G M (strict convexity). -i 
The conditions of Definition|4]resemble the conditions of Definition|2l particularized to horizontal 
tangent vectors 5x = T x Ai \ V X = H x \ {0}. The induced metric © is only a pseudo-distance 
on Ai since two states x a , x b £ Ai may satisfy d(x a , x b ) = despite x a ^ x b . In fact, every 
piecewise differentiable curve 7: I — >Ai that satisfies 7(5) G V 7 ( s ), for almost every s £ I, also 
satisfies that L V(7(s), j(s))pds = 0. The pseudo-distance d measures the "distance" between 
two given points x a and x b by considering only the horizontal component of curves 7 : / — > Ai 
connecting x a and x b , that is, the component 7ft(s) of 7(5) = 7/i( s ) + iv( s ) where 7h(s) G % 7 ( s ) 
and 7„(s) G V 7 ( s ), for each s £ I. 

We can now provide the reformulation of Theorem[T]for horizontal Finsler-Lyapunov functions. 

Theorem 3: Consider the system (OQ) on a smooth manifold Ai with / of class C 1 , a vertical 
distribution V x (|42l) . and a horizontal distribution H x (|43T ). Let C C M be a connected and 
forward invariant set, C Ai a set of isolated points, and a a function in M> — )■ R>o- 
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Given a candidate horizontal Finsler-Lyapunov function V for (OQ) on H x , suppose that © 
holds for each t > t , each x G C and each 5x G 7^. Then, the solutions to (QQ) 

(i) do not expand the pseudo-distance d © on C if a(s) = for each s > 0: for any xi, x 2 G C 

and any t > t , d(ip to (t,x 1 ),ip to (t,x 2 )) < j(d(x 1 ,x 2 )), where 7(5) > s; 

(ii) asymptotically contract the pseudo distance d on C if a is a /C function: (i) holds and for 
any x u x 2 G C, lim d{ip tQ {t,Xi),ip to {t,x 2 )) = 0; and 

t— >oo 

(iii) exponential contract the pseudo distance d on C if a(s) = As > for each s > 0: there exists 

K > 1 s.t., for any 07, x 2 G C and any t > t , d(ip to (t, ^t (^ ^2)) < Ke^ x ^~ to ^d(xi, x 2 ). 

j 

Proof: [Sketch] Following the proof of Theorem \T\ for any given curve 7 : 7 — )■ C, let 
^t ('5 7( s )) be the solution to (Q~|) from the initial condition 7(5). We have to prove that the length 
of the curve V*o(*>7(0) S iven b y // ^(^to(*> 7(*))> D ^t {t, 7(*))[0> 7(s)])*ds is decreasing (not 
increasing) as t — > 00. Define x(t, s) := ^(t, 7(5)), and &r(t, s) := Di/) to (t, j(s))[0, S'y(s)]; 
and consider the decomposition of &c(t,s) into Sx(t,s) = 5xh(t,s) + Sx v (t,s), respectively 
horizontal 5xh(t,s) G H x (t, s ) and vertical Sx v (t, s) G V x (t )S ) components. 

From the assumption (ii) on V, V(x(t, s), 5x(t, s)) = V(x(t, s), 5x h (t, s)), thus 

-f-Vfxft, s),6x(t,s)) = ^-V(x(t,s),8x h (t,s)) for each t > t , s G I. (44) 
at at 

From the assumption (ii) on V, since © holds for each t > t , x G C and 5x G 7^, by repeating 

the argument in the proof of Theorem [H for each t > t and s £ I we get 

s), 6x h {t, s)) < -a{V{x{t, s), Sx h {t, s))) = -a(F(:r(t, s), <fcr(t, s))). (45) 

at 

Combining (O and ((45]) we get the inequality ^K(x(t, s), 5x(t, s)) < —a(V(x(t, s), Sx(t, s))) 
that coincides with (fl"5l) . The result of the theorem is then established by following the proof of 
Theorem OQ from CE2>. ■ 
Remark 4: The formulation of the LaSalle-like relaxations of Theorem [2] and Remark [3] in 
Section IVIIl immediately extend to horizontal Finsler-Lyapunov functions. Following Remark [2l 
the regularity assumption (i) in Definition |4] can be relaxed to functions V that are piecewise 
continuously differentiable and locally Lipschitz. In such a case, the goal is to show that the 
inequality (fT8l) holds. This is guaranteed, for example, if (1431) holds for almost every t and s. j 
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B. Horizontal contraction and contraction on quotient manifolds 

The notion of horizontal space is classical in the theory of quotient manifolds. Let M. be 
a given manifold and let Ai\ ~ be the quotient manifold of M. induced by the equivalence 
relation ~G M. x M.. Given x E M, we denote by [x] GM\~ the class of equivalence to x. 
Suppose that the system x = f(t,x) in CD is a representation on M. of a system on .M\~ in 
the following sense: for every t > 0, every x , and every z E [x ], the solution z ) to (OQ) 
satisfies y>t (t, £o) [Vto(^) x o)] f° r eacn £ > V I n sucn a case we cai l ^ = f{tj x ) a quotient 
system on .M\ ~. The equivalence relation ~ usually describes the symmetries on the system 
dynamics on A4, which implicitly characterize the quotient dynamics. Every solution Lf to (-,z ) 
of (OQ) from z E [x ] E A^\~ is a (lifted) representation of a unique solution [ip to (-,xo)] on the 
quotient manifold. 

The vertical space V x at x is defined as the tangent space to the fiber through x. In this 
way, any tangent vector 5[x] to TmM\ ~ has a unique representation in the horizontal space 
T-L x , called the horizontal lift [1J. The particular selection of the vertical distribution guarantees 
that the horizontal Finsler-Lyapunov function V on T-L x is zero for each 5x E V x , that is, V 
is a Finsler-Lyapunov function for the quotient system on M.\ ~. The induced metric d © 
is a pseudo-distance on A4, but utf zie \ Xl ] )Z2e \ X2 \d(zi,Z2) is a distance on A^\~ for any given 
[xi], [x 2 ] E AA\~. Therefore, establishing asymptotic contraction of the pseudo-distance d for 
© on M. is equivalent to establish incremental asymptotic stability of the quotient system on 
.M\~, implicitly represented by (OQ) on Ai. 

Suppose that Theorem |3] holds for a given quotient system (OQ). Then, by considering the lifted 
solutions of (OQ) to A / l\~, the system (OQ) is (i) incrementally stable on C if a(s) = for each 
s > 0; (ii) incrementally asymptotically stable on C if a is a K, function; and (in) incrementally 
exponential stable on C if a(s) = As > for each s > 0. In this sense, horizontal contraction 
in the total space is a convenient way to study contraction on quotient systems. 

Quotient systems are encountered in many applications including tracking, coordination, and 
synchronization. The potential of horizontal contraction in such applications is illustrated by two 
popular examples. 

Example 5: [Consensus] 
We consider consensus algorithms of the form 

x = A{t)x (46) 
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where and, for each t > 0, A(t) has nonnegative off-diagonal elements and row sums zero 

(we assume that A(t) is continuously differentiable). These Metzler matrices |fT9l are typically 
used to model the graph topology of network problems. Indeed, the 5-graph ofA{t) has an arch 
from the node i to the node j, i ^ j, if > S > 0. 

Given 1 := [i ... i] T , the row sums equal to zero guarantee that A(t)l = for each t > 0. 
Indeed, al is a consensus state of the network for every a G R. Because of this symmetry, (|46|) 
represents a quotient system on the quotient manifold R n \ ~ constructed from the equivalence 
x ~ y iff x — y = al, for some a > 0. In fact, if x ~ y then = A{t)y for each t > 0. The 

elements of R n \ ~ are [x] := {x + al\a e R}, the vertical space is given by V x := Span({l}), 
and the horizontal space can be taken as 1-L X := {Sx G R n | l T Sx = 0} = V^r. 

To show contraction, consider the displacements dynamics from (l46l) given by 5x = A{t)8x, 
and the horizontal Finsler-Lyapunov function 

V(x, 5x) := max 5xi — min Sxi, (47) 

i i 

that coincides with the classical consensus function adopted in |fT9l , ll36l lifted to the tangent 
space. Note that V satisfies every condition of Definition @] but continuous differentiability. In 
particular, V is positive and homogeneous for every Sx G T-L x . For Sx G T x R n , V(x, Sx) = 
V(x, Sxh) with Sxh horizontal component of Sx, since V(x, Sxh + al) = V(x, Sxh) for each 
a G R. 

Following Remark HI the lack of differentiability is not an issue. In fact, from |fT9l Section 
3.3], for any initial condition x G R n and any initial tangent vector 5x G T x R n = R n , V is 
non-increasing along the solution (p to (-,x ) to (|46|) . namely V(ip to (t, x), D(p tQ (t, x )[0, Sx Q ]) < 
V(x ,Sx ) for each t > t . This inequality is the result of the combination of [19, Section 
3.3], showing that maxzj — minzj is non-increasing for z = A{t)z, and of the fact that 

i i 

the evolution Dip to (t, x )[0, Sx ] of Sx along the solution cp to (-,Xo) is also a solution to the 
differential equation Sx = A(t)Sx (as shown in (fT4l)). 

By the same argument, exponential decreasing of V is achieved under additional conditions on 
uniform connectivity on the adjacency matrix A(t). Following [fT9l Theorem 1], define A*(t) : = 
f* +T A(r)dT and suppose that there exist k G {1, . . . , n}, S > 0, and T > such that, for every 
t > to and every j G {1, . . . ,n} \ {k}, there is a path from the node k to the node j of the 
5-graph of A*{t). Then V decreases exponentially along the solutions to (|46|) . By integration, the 
quotient system defined implicitly by (|46|) is incrementally exponentially stable. As a corollary, 
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every solution to the quotient system converges to the steady-state solution [0], that is, every 
solution to (1461) exponentially converges to consensus. 

The reader will notice that the incremental exponential stability of (l46l) is a straightforward 
consequence of the exponential stability results of lfT9l , through the lifting to the tangent space 
of the (non-quadratic) Lyapunov function used in |fT9~l . In this sense, the differential framework 
captures the equivalence on linear systems between stability and incremental stability. j 

Example 6: [Phase Synchronization] 
Consider the interconnection of n agents 9 k = u k , 9 k G S 1 (phase), given by 



0u 



n 

-Y 



sin 



0i 



(48) 



Using Sj k := sm(9j — 9 k ), c jk := cos(9j — 9 k ), 1 := [i ... i] T , the aggregate state 9 := [o- L ... e n ] T , 
and the displacement vector 59 := [<S0i ... se n ] T , (|48l and the related displacement dynamics can 
be written as follows. 



1 

n 



r o S2i 

S12 o 



Snl 
Sn2 



59 



n 



:S(fl) 



- J2 C J1 c 21 

Cl2 - J2 c j2 

3^2 



Cnl C„2 



Cn2 



c i" 



59. 



(49) 



= :C(0) 



(1491 ) is a quotient system based on the equivalence 9 ~ iff there exists a G R such that 
6>-6> = la. In fact, S(6>) = S(6>+al), which fixes the class of equivalence [9] = {9+al \ a G R}, 
and the vertical space Vq := Span{l}. As in the previous example we consider T-Lg := Vq = 
{59 G R" 1 1 T 59 = 0}. 

It is well-known that the open set O C E> n , given by phase vectors 9 such that \9j — 9 k \ < 3 for 
each j,k G {1, . . . , n}, is forward invariant. Therefore, for any given initial condition 9 G O, 
the solution ip(-,9o) to (|49|) guarantees that C(ijj(t, 9q)) is a Metzler matrix for each t > 0. It 
follows that (gU) contracts on the lifted Hilbert metric V(0, 59) := max 5^ - min 59, (|47l) . 

i i 

Moreover, since C(#) is symmetric, (l49l) contracts the horizontal Finsler-Lyapunov function 
^(5,50) := 5^ T [/ n - ^]59, as shown in E3 Proposition 1]. 

Paralleling Example [3l almost global contraction can be established by considering the hori- 
zontal Finsler-Lyapunov function given by the non-flat metric 

1 



V(9,59) :-- 



P 



2q 



59 1 U59, 



qeN, 



(50) 
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where II :- 



J _ ill 



(note that U89 = for 89 E V$), and p is the magnitude of the centroid 
p e %<t> ■= - e' l9k . Following Il29"1 , p E [0, 1] is a measure of synchrony of the phase variables, 

n k=l 

since p is 1 when all phases coincide, while p is when the phases are balanced, p is also 

n 

nondecreasing, since p = £ sin(9 k — 4>) 2 . In particular, p = for p = (balanced phases) 

n fc=i 

n 

or for ^ sin(# fc — <p) 2 = 0, which occurs on isolated critical points given by n — m phases 

k=l 

synchronized at <p + 2jrc and m phases synchronized at <p + ti + 2jir, for j 6 N and < m < | 
Synchronization is achieved for m = 0, the other critical points are saddle points (for an extended 
analysis see [29, Section III]). 

Using V to denote the left-hand side of ©, we get 

v = jr q te T (- 2 -± Esin(^-0) 2 n + nc(e) + c(^n^^ 



(51) 



For each G S", V" = for 5# e V e . V" is negative for 9 e O and 5^ e "Hg. For G §" \ 
and 56* G Hg, g can be suitably chosen to balance the presence of positive eigenvalues in C(6>). 
In fact, given any compact and forward invariant set C C §" that does not contain any balanced 

n 

phase (p = 0) or saddle point sin(# fc — 0) 2 = 0), there exists a sufficiently small e > such 

k=l 

n 

that sin(6'fc — cp) 2 > e and p > for every 9 E C. Thus, contraction on C is established by 

k=l 

picking q > -. j 



C. Forward contraction 

The use of horizontal contraction is not restricted to quotient systems. We briefly discuss 
in this section the concept of forward contraction of x = f(x), that we define as horizontal 
contraction for the particular case 

H x := Span({/(x)}), for each x E M. (52) 

By definition, forward contraction captures the property that for every solution <p(-,xo) to x = 
f(x), x E M, and every T > 0, the points <p(t + T, x ) and <p(t, x ) converge to each other as 
t — > oo. This property has strong implications for the limit set of x = f(x), as illustrated by the 
following proposition. Restricting the analysis to time-invariant systems x = f(x) for simplicity, 
we propose two novel results on attractor analysis by exploiting forward contraction. 
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Proposition 2: [Bendixon's like criterion] Consider the system i = /(x) on a smooth manifold 
Ai with / of class C 1 , and a forward invariant set C C At. Given a K, function a and a candidate 
horizontal Finsler-Lyapunov function on % x in (l52l) . suppose that Theorem [3] holds for ± = f(x). 
Then, no solution of x = f(x) in C is a periodic orbit. j 
Proof: Suppose that from x E C, the solution </?(•, x ) is a periodic orbit T. Then, from 
the definition of % x and the continuity of V, there exist m > such that m < V(x, /(x)) for 
each x G T (r is a compact set). From ©, the definition % x , and the fact that a is a function 
of class /C, there exists a class ICC function (3 such that m < lim V(ip(t, x Q ), f(ip(t,x ))) < 

t— >oo 

lim l3(V(i>(0, x ), f(ip(Q, x ))),t) = 0. A contradiction. ■ 

t— >oo 

Forward contraction makes contact to a vast body of theory, primarily motivated by the 
Jacobian conjecture [6J. Conditions to establish the absence of periodic orbits are proposed in 
ll32l (see e.g. Theorem 7) and Il20l . and are based on specific matrix measures. The connection 
to Theorem Q] can be established along the lines of Section |VIJ These conditions are generalized 
in ifTTl . which connects the absence of periodic orbits to the contraction of a suitably defined 
functional S in the manifold tangent bundle, as shown in Q/71 Sections 2 and 3]. In a similar way, 
Proposition |2] relates the absence of periodic orbits to the contraction of a horizontal Finsler- 
Lyapunov function V on % x = Span{/(x)}. 

Under the assumption of boundedness of the solutions to x = f(x), the absence of periodic 
orbit induced by the contraction argument is exploited in the next proposition to guarantees that 
a given set A is asymptotically attractive. 

Proposition 3: [Asymptotic attractor on C] Consider the system x = f(x) on a smooth 
manifold M. with / of class C 1 , a forward invariant set C C M., and a forward invariant set 
(attractor) A C C. Given a /C function a and a candidate horizontal Finsler-Lyapunov function 
on H x in (|52|) . suppose that Theorem [3] holds for x = f(x), with the relaxed condition that © 
holds for each x E C \ A, and each 5x E T-L x . If 

• A contains every equilibrium point = f(x), x E C; 

• for every initial time t and every initial condition x E C, there exists a bounded set 
Mx ^ M. such that ip{t, x ) E U XQ for each t > 0, 

then A is a forward invariant set, and for every initial condition x E C, and every neighborhood 
U D A, there exists T^, ^ > such that ip(t,x ) E U for each t > T{x ,u)- j 
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Proof: Since ip(t,x Q ) belongs to the bounded set U xo for each t > 0, by [fT51 Lemma 4.1] 
it converges to its w-limit set, given by the compact and forward invariant set lu + (x ) := {x G 
M. | x = lim i/j(t n ,x ) where t n G M>o — >■ oo as n — > oo}. Note that if lim i/j(t,x ) = x* G C 

n—>oo t—too 

then, by hypothesis, x* belongs to A C U. Therefore uj + {xq) \ A does not contains equilibria. 
We prove by contradiction that tu + (x ) C A. 

Suppose that cu + (x ) fl A = 0. By compactness of w + (a;o), the definition of H x , and the 
continuity of V, there exist m > such that m < V(x,f(x)) for each x G w + (x ). Con- 
sider the solution ip(-,x) whose initial condition x G co + (x ). From ©, the definition T^., 
and the fact that a is a function of class /C, there exists a class XX function such that 
m < lim ^(^(t, x), f{ip(t, x))) < lim /3(y(-0(O, x), f{ip{0, x))),t) = 0. A contradiction. 

Suppose that cu + (x )r}A ^ and cu + (x ) <£. A. By the same argument used above, there exists 
a sequence of t k G M>o such that t k — > oo as k — > oo such that V(-ip(t k , x), f(-ip(t k , x))) > m > 
but lim V(ip(t k ,x),f(^(t k ,x))) < lim /3(V(i/j(0, x), f(i/j(0, x))), t k ) = 0. A contradiction. ■ 

IX. Conclusions 

The paper introduces a differential Lyapunov framework for the analysis of incremental 
stability, a property of interest in a number applications of nonlinear systems theory. Our main 
result extends the classical Lyapunov theorem from stability to incremental stability by lifting the 
Lyapunov function in the tangent bundle. In addition to classical Lyapunov conditions, Finsler- 
Lyapunov functions induce a Finsler differentiable structure in the state space. Through inte- 
gration along curves, the construction of a Finsler-Lyapunov function, a local object, implicitly 
provides the construction of a decreasing distance between solutions, a global object. 

The study of global distances through local metrics is the essence of Finsler geometry, a 
generalization of Riemannian geometry. Several examples and applications in the paper suggest 
that the Finsler differentiable structure is indeed the natural framework for contraction analysis, 
unifying in a natural way earlier contributions restricted either to a Riemannian framework lfl8ll . 
or to matrix measures of contraction E71 . Il33l . In the same way, the formulation of the 
results on differentiable manifolds rather than in Euclidean spaces is not for the mere sake 
of generality but motivated by the fact that global incrementally stability questions arising in 
applications involve nonlinear spaces as a rule rather than as an exception. 

A central motivation to bridge Lyapunov theory and contraction analysis is to provide contrac- 
tion analysis with the whole set of system-theoretic tools derived from Lyapunov theory. The 
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present paper only illustrates this program with LaSalle's Invariance principle but we expect 
many further generalizations of Lyapunov theory to carry out in the proposed framework. This 
includes the use of asymptotic methods such as averaging theory or singular perturbation theory 
(see e.g. the result ), and, most importantly, the use of contraction analysis for the study of 
open and interconnected systems. The original motivation for the present paper was to develop 
a differential framework for incremental dissipativity flU, [fT2l . 041 . which will be the topic of 
a separate paper. 

Although a straightforward extension of contraction, the concept of horizontal contraction 
introduced in this paper illustrates the potential of contraction analysis in areas only partially 
explored to date. Primarily, it provides the natural differential geometric framework to study 
contraction in systems with symmetries, disregarding variations in the symmetry directions where 
no contraction is expected. Problems such as synchronization, coordination, observer design, 
and tracking all involve a notion of horizontal contraction rather than contraction. The notion of 
forward contraction, which corresponds to the particular case of selecting the vector field to span 
the horizontal distribution, connects the proposed framework to an entirely distinct theory which 
seeks to characterize asymptotic behaviors by Bendixon type of criteria, excluding periodic orbits 
or forcing convergence to equilibrium sets [17J. Finally we suspect that the density approach 
to stability analysis proposed by Rantzer E51 is closely connected by duality to the proposed 
differential framework, and this will also be the object of future studies. 

Overall, we anticipate a number of interesting developments beyond the basic theory presented 
in this paper and we hope that the proposed differential framework will facilitate further bridges 
between differential geometry and Lyapunov theory, a continuing source of inspiration for 
nonlinear control. 
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